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Abstract. The times of maximum brightness collected in the GEOS RR Lyr database al- 
lowed us to trace the period variations of a sample of 123 galactic RRab variables. These 
data span a time baseline exceeding 100 years. Clear evidence of period increases or de- 
creases at constant rates has been found, suggesting evolutionary effects. The observed rates 
are slightly larger than those predicted by theoretical models; moreover, there is an unex- 
pected large percentage of RRab stars showing a period decrease. The new possibilities 
offered by the use of robotic telecopes (TAROTs, REM) and of data from satellite (CoRoT) 
are expected to speed up the project to measure stellar evolution in real time. 
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1. Introduction 

Horizontal branch stars cross the instability 
strip in different stages of their evolution; in 
these phases they become RR Lyr stars. The 
crossing can take place in both directions; as 
a consequence, the periods will be either in- 
creasing if the stars evolve from blue to red 
or decreasing if they evolve from red to blue. 
Despite its importance as a test for the stellar 
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evolution theory, the observed rate of the pe- 
riod changes is still an unknown quantity. 

The amateur/professional association 
GEOS (Groupe Europeen d' Observations 
Stellaires) has built a database aiming to put 
together all possible RR Lyr light maximum 
times published in the literature. The publi- 
cations from the end of the XIX th century up 
to today were scanned for this purpose and 
recent observations from amateur astronomers 
of the European groups GEOS and BAV 
(Bundesdeutsche Arbeitsgemeinschaft fur 
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Fig. 1. Histogram of the ratios between the modulus 
of the a coefficient and its formal error. 

Veranderliche Sterne) fed the database in a 
continuous way. To date, it contains about 
50000 maximum times from more than 3000 
RR Lyr stars. The GEOS database is freely 
accessible on the internet at the address 
http://dbrr.ast.obs-mip.fr/. Therefore, this 
database is the adequate tool to search for 
secular period variations in the galactic RR Lyr 
stars, since most of them have been studied 
for tens of years, several of them since the 
end of the XIX th century. The first result of 
the project, Stellar e volution in real t i me, hav e 
been presented by iLe Borgne et al] (l2007al) . 
It is noteworthy that its outlines have been 
sketched during several GEOS meetings, 
where the different knowledge of amateur 
and professional astronomers found a very 
profitable synthesis. 

2. The methodological approach: 
some insights 

To calculate the rates of the period changes, the 
times of maximum brightness have been fitted 
by means of the linear elements 

C = T + PE (1) 

where P is the period (in days) and E is the cy- 
cle counter. In several cases the linear fit was 
not able to reproduce the observed times of 
maximum brightness and a parabolic trend ap- 
peared in the O-C (observed minus calculated 



values) plots. Therefore, the maxima were fit- 
ted with a second-order polynomial 

C = T + P ■ E + a ■ E 2 (2) 

where Pq is the period at epoch Tq. The value 
of the a coefficient is related to the rate of the 
period change by means of 

dP/dt = 2 • a/(P) (3) 

where dP/dt is measured in [d/d] and (P) 
is the average period, i.e., the ratio between 
the elapsed time and the elapsed cycles. 
Theoretical models often refer to the rate of pe- 
riod changes using the /3 parameter, expressed 
in [d Myr~'], since it directly indicates how 
much the period change in 10 6 years. 

Though the times of maximum brightness 
have been collected by using very different 
techniques (visual, photographic, photolectric, 
CCD), the very long time baseline allowed us 
to trace the period variation in a clear way. 
We extracted from the GEOS database the 
RRab stars for which the times of maximum 
brightness span more than 50 years. Fifty-four 
stars show a constant period, while 27 show 
an increasing period and 21a decreasing one. 
The remaining 21 RRab variables show erratic 
changes. Figure Q] shows the histogram of the 
ratios between the modulus of the a coefficient 
and its formal error for the 48 stars showing 
a linearly variable period. The ratios are al- 
ways greater than 5.0 and are greater than 10 in 
all cases except four. Therefore, the parabolic 
trends are very well defined and they bear wit- 
ness of the reliability of the results. 

Figure [2] illustrates some examples of the 
O-C behaviours. Among the stars showing a 
period increase (top row), BN Aqr is one of 
the most evident cases: the parabolic trend left 
by the linear elements is very well defined 
and the two branches are symmetrical. No im- 
provement has been obtained by introducing a 
third order term; as a matter of fact, the resid- 
ual O-C values are very flat around the zero 
value. The O-C plot from the linear fit of the X 
Ari maxima is slightly different: the ascending 
branch seems to be steeper than the descend- 
ing one. As a consequence, the residual O-C 
values show an apparent oscillation. We care- 
fully investigated the nature of this oscillation, 
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Fig. 2. For each star the upper panel shows the O-C values as obtained from linear elements, the lower 
panel from a parabolic fit. Measure units are elapsed cycles (E/10000) on the x-axis and days on the y-axis. 



visible also in other cases. The related period- 
icities are often about half the time coverage 
and their amplitudes are quite small. Therefore, 
we can infer that the parabolic fit is not perfect 
and that some systematic residuals were left. 
The maxima of RS Boo shows this effect in a 
more pronounced way. In some cases we im- 
proved the O-C fitting by calculating a third- 
order least-squares fit, but in three cases only 
we obtained a significant third-order term. We 
note that no appreciable change is derived for 
the value of the a coefficient when performing 
a third-order fit. The case of RS Boo is still 
more intriguing, since the residual O-C val- 
ues show a scatter on another timescale, shorter 
than the two previous ones. As a matter of fact, 
when performing the frequency analysis of the 
residual O-C values (i.e., the prewhitened data 
after the third-order fit) we were able to de- 
tect a peak corresponding to at 534 d (Fig. [3]», 
practically coincid ent with the value given by 
lOosterhooild 19491) for the Blazhko period, i.e., 
533 d. The connection between period changes 



and the Blazhko effect is still unclear and de- 
serves more attention in the future. 

Among the stars showing a period decrease 
(bottom row), V964 Ori shows the clearest pat- 
tern, in spite of the small number of points. 
Asymmetrical O-C plots are also common, 
perhaps more common than in the case of 
period increases, and VX Her constitutes a 
good example. As in the case of RS Boo, also 
AH Cam shows a Blazhko effect, but on a 
very s hort timescale (only 1 1.18 d, Smith et al. 
1994): it is notewothy that also in this case the 
times of maximum brightness are able to deter- 
mine the correct Blazh ko period (see Fig. 6 in 
iLe Borgne et~alll2007ah . 

3. Scientific conclusions 

The most remarkable result is the detection of a 
large number of O-C plots showing a parabolic 
pattern when linear elements are used. The 
cause of this pattern is a regular period varia- 
tion, having the same order of magnitude of the 
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Fig. 3. RS Boo: power spectrum of the residual 
O-C values after having fitted the maxima with 
a third-order polynomial. The highest peak is at 
0.0019 d -1 and it corresponds to the known period 
of the Blazhko effect. 

expected evolutionary changes. On the other 
hand, in 54 cases (i.e., a number of cases simi- 
lar to that of RR Lyr stars showing period vari- 
ations) we did not detect a significant period 
variation and linear elements satisfy the avail- 
able times of maximum brightness. The upper 
panel in Fig. H] shows the distribution of the ft 
values in the sample of the galactic RRab stars. 
The central peak is given by the stars showing 
a constant period, but the number of stars both 
in the negative and in the positive planes are 
relevant. 

Indeed, the comparison between RRab 
stars showing negative and positive rates yields 
us the most important scientific feedbacks. 
They can be evaluated in an easier way by 
removing the central peak (lower panel in 
Fig. 0}. It is quite evident that the blueward 
evolution (decreasing periods, negative fj val- 
ues; 21 cases) is as common as the redward 
one (increasing periods, positive fj values; 27 
cases). B y analyzing and c ompari ng these two 
samples, iLe Borgne et al.l (l2007al) could stress 
some well-established observational facts: 

1. RRab stars experiencing blueward evolu- 
tion are quite common, only slightly less 
than RRab stars experiencing redward evo- 
lution. The period range covered by the two 
groups are very similar and the average and 
median periods are nearly coincident; 

2. the absolute values of period changes are 
larger than expected. The t heoretical e volu- 
tionary rates calculated by lLed (Il99ll) pre- 
dict fj values smaller than 0.30 for almost 
all stars. fj values greater than 0.30 are re- 
ported in exceptional cases. Since in our 
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Fig. 4. Histograms of the distributions of the /? val- 
ues of the whole sample (upper panel) and restricted 
to the stars shown an appreciable linearly variable 
period (lower panel). The observed /? values are 
slightly larger than those predicted by the evolution- 
ary models. 

sample we have 18 stars with 0.00 < fj < 
0.30 and 9 stars with fj >0.30, the 2:1 ratio 
seems too small to be caused by the excep- 
tional conditions. 
3. Blazhko effect is often superimposed on 
secular changes, but the monotonic trend 
due to evolutionary variations still remains 
visible; 

In particular, the rates of the period changes are 
difficult to associate with canonical horizontal- 
branch evolution. The large observed fj values 
could be ascribed to a sort of "noise" super- 
imposed on the canonical evolutionary tracks, 
irrespective of where the star is in the course 
of its horizontal-branch evolution. However, 
this is still surprising since the main sources 
of this "noise" had been invoked to act toward 
the exhaus tion of the helium in the c ore, but 
not before dSweigart & Renziniiri979t) . At this 
purpose, we also note that there is a numer- 
ous group of RRab stars showing erratic O-C 
behaviours, with abrupt or continuous period 
variations, never monotonic. The regular vari- 
ations caused by light-time effect (and hence 
duplicity), often invoked to explain large O- 



Poretti et al.: Stellar evolution in real time 



5 



C excursions, are not convincingly observed 
in any case (except maybe for RZ Cet) of our 
sample. 

4. The future 

As a general conclusion about the comparison 
between our observational results and the the- 
oretical predictions, we claim that there is a 
very powerful feedback between the two ap- 
proaches. In particular, theoretical investiga- 
tions should take into account that many RR 
Lyr stars are showing a blueward evolution. 
The theoretical models should also match the 
observed f3 values in a more satisfactory way, 
as these seem to be higher than expected, both 
for red ward and blueward evolutions. 

To rapidly progress in the project Stellar 
evolution in real time we need to accumu- 
late much more accurate times of maximum 
brigthness. In such a case the parabolic trend 
will be evidenced over a time baseline much 
shorter than the 110 yr interval used for this 
first approach. The robotic telescopes TAROTs 
(Telescope a Action Rapide pour les Objets 
Transitoires or Ra pid Action Telescop e for 
Transient Objects; Brin ger et al.1 1 19991) . op- 
erating at La Silla (Chile) and at Calern 
Observatory (France), are currently used to 
monitor RRab stars brighter than V=13.0. 
They are providing very homogenous list of 
CCD maxima (e.g., iLe Borgne et al J l2007bl) 
characterized by small errorbars (+0.003 d on 
average). 

Another critical item is to investigate the 
interplay between Blazhko periods and evolu- 
tionary changes, which can be evidenced only 
by means of dense lists of maxima over one or 
two decades. We also have to note that the driv- 
ing mechanism of the Blazhko phenomenon 
needs some clarification, since the two com- 
petitive explanations, the magnetic and the res- 
onance models, are both still plausible. To pur- 
sue both goals in a specific way, we started 
a survey of few Blazhko RRab stars by us- 
ing the robo t ic RE M (Rapid Eye Movement, 
ICovino et al.1 120031) telescope located at La 
Silla, Chile. 

Moreover, a substantial improvement in the 
study of the light curves of RR Lyr variables 



could be obtained by a continuous monitor- 
ing, without the seasonal or the night-day al- 
ternances. Such a possibility was a dream un- 
til a few years ago. After the successful launch 
on December 27th, 2006, the satellite CoRoT 
(COnvection, ROtation and planetary Transits) 
is now transforming the dream into reality. 
Staring on the same field for 150 d, CoRoT 
will monitor 10 bright stars in the asteroseis- 
mologic channel and 12000-15000 in the exo- 
planetary one. In the latter channel, the moni- 
toring of some classical variables can be pro- 
posed in the framework of the 'Additional 
Programs" and a specific one on RR Lyr vari- 
ables has been approved (PL: M. Chadid, Nice 
University; some of the authors are Cols). In 
such a context, the monitoring of three RR 
Lyr variables (VI 127 Aql, V1220 Aql, and 
NSVS 12568727) has been proposed for the 
two first CoRoT runs with the goal to investi- 
gate very subtle cycle-to-cycle variations (and 
hence large or small Blazhko effects, if any). 
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